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1.  Introduction 


The  general  goal  of  this  research  program  was  to  study  the  various  physical 
mechanisms  involved  in  the  onset  of  leading-edge  separation  and  vonex  formation  when 
airfoils  execute  highly  unsteady  motions.  Unsteady  separation  and  dynamic  stall  are 
inherently  complex  processes  which  are  influenced  by  many  interrelated  flow  parameters. 
Instead  of  a  parametric  study,  our  investigation  attempted  to  isolate  the  basic  fluid 
mechanical  phenomena  for  the  specific  case  of  incompressible  flow  around  a  NACA-0012 
airfoil  pitching  (about  the  1/4-chord  axis)  to  large  angles  of  attack  (typically  60  degrees). 
Both  constant  pitch  rate  (i.e.  ramp-type)  motions,  and  variable  pitch  rate  motions  were 
studied. 

Our  work  initially  started  as  purely  an  experimental  investigation  and  quickly  lead 
to  a  coordinated  experimental  and  computational  effort.  The  computational  work  -  a 
collaboration  with  Dr.  Miguel  Visbal  (WPAFB)  -  allowed  us  to  obtain  information  on  the 
flow  variables  that  were  difficult  to  measure  (such  as  vorticity  and  reversed  flow  evolution 
near  the  surface)  for  our  experiments.  This  close  collaboration  between  experiment  and 
computation  has  proved  to  be  an  ideal  approach  in  the  study  of  such  complicated 
phenomena  as  leading  edge  separation  and  dynamic  stall. 

This  final  report  provides  a  brief  sunvnary  of  the  highlights  of  our  results  to  date. 
Most  of  these  results  have  been  described  in  much  greater  detail  in  previous  annual  reports 
[1-4]  and  conference  presentations/papers  [5-7].  Some  of  the  results  have  recently  appeared 
in  a  journal  publication  (see  Appendix  A);  one  paper  is  currently  under  review  (see 
Appendix  B).  Two  other  manuscripts  are  being  prepared  for  publication.  This  report  will 
conclude  by  discussing  some  of  the  outstanding  issues  which  need  further  research. 


1.  Summary  of  Results  for  Constant  Pitch  Rate 

The  experiments  and  computations  for  this  case  involved  very  rapid  acceleration  to 
a  constant  pitch  rate  trajectory.  The  chord  Reynolds  number  was  of  order  10^  (8,000  ^  Re^ 
^  15,000)  and  non-dimensional  pitch  rate  varied  between  0.015  ^  fi*  <  0.4.  The  main 
findings  were: 

■  The  onset  of  leading  edge  separation  can  be  determined  from  flow  visualization 
experiments  by  marking  the  appearance  of  a  kink  or  a  "bulge"  near  the  leading 
edge.  Computations  show  that  this  feature  correlates  very  well  with  the  appearance 
of  a  similar  feature  in  the  vorticity  field.  Visual  and  numerical  predictors  of  the 
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separation  angle,  determined  from  the  computations,  agree  well  with  our 
experimental  results. 


■  The  shear  layer  that  forms  between  the  region  of  reversed  flow  near  the  airfoil 
surface  and  the  outer  inviscid  stream  appears  to  behave  like  a  "free"  shear  layer. 
It  is  subject  to  Kelvin-Helmholtz  instability  and  exhibits  many  of  the  vortex 
interactions  (such  as  pairing)  observed  in  free  shear  layers.  Since  a  counter  flow 
situation  exists,  this  shear  layer  may  exhibit  absolute  (as  opposed  to  convective) 
instability  characteristics  under  some  conditions. 

■  At  low  pitch  rates,  onset  of  leading  edge  separation  occurs  after  the  reverse  flow 
originating  from  the  trailing  edge  has  reached  the  leading  edge  (see  Figure  1). 

■  At  high  pitch  rates,  flow  separation  is  leading-edge  dominated;  a  region  of  reverse 
flow  develops  near  the  leading  edge  independently  from  that  originating  from  the 
trailing  edge  (see  Figure  2). 

■  Various  estimators  for  determining  the  onset  of  leading  edge  separation  were 
compared.  These  estimators  were  based  on  the  behavior  of  the  vorticity  Held,  the 
reverse-flow  evolution  near  the  surface,  surface  pressure  and  pressure  gradient  fields. 
All  estimators  showed  the  same  trend  for  separation  angle  versus  non-dimensional 
pitch  rate.  The  estimator  based  on  the  evolution  of  reverse-flow  appears  to  hold  the 
best  promise  for  flow  control  applications. 


2.  Summary  of  Results  for  Non-Constant  Pitch  Rate 

Studies  of  non-constant  pitch  rate  included  two  types  of  motions.  Brst  were  airfoil 
trajectories  which  deviated  from  the  ideal  ramp  due  to  rinite  acceleration  and  deceleration 
periods.  Second  were  variable  pitch  trajectories  constructed  from  three  constant  pitch  rate 
segments  with  non-dimensional  pitch  rates  of  0.5f2*,  and  212*. 

(a)  Effect  of  Initial  Acceleration 

■  The  elapsed  time  for  onset  of  leading  edge  separation  is  affected  by  the  initial 
acceleration. 
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■  The  angle  of  attack  where  separation  is  observed  is  not  sensitive  to  initial 
acceleration  for  the  relatively  high  pitch  rate  motions  studied. 

■  Many  of  the  details  of  dynamic  stall  vortex  formation  and  its  interactions  (e.g. 
convection  speed  11^)  appear  to  be  also  unaffected  by  the  initial  acceleration. 

■  The  initial  acceleration  was  observed  to  affect  the  integrated  loads,  surface  pressure 
distribution,  and  the  evolution  of  the  reverse  flow  regions  only  during  the 
acceleration  period  and  for  a  relatively  short  time  (t  s  tU^C  ^  0.25)  afterwards. 
After  that  time,  all  of  these  quantities  depend  only  on  the  instantaneous  angle  of 
attack  for  a  given  pitch  rate. 


(b)  Variable  Pitch  Rate  Results 

■  The  instantaneous  pitch  rate  a,  ^  at  the  static  stall  angle,  believed  to  be  important 
for  the  case  of  oscillating  airfoils,  was  found  not  to  be  a  controlling  parameter  in 
general. 

■  Our  working  hypothesis  implies  that  if  the  airfoil  motion  time  scale  is  fast,  onset 
of  leading  edge  separation  is  expected  to  be  based  on  the  current  pitch  rate  (and 
angle  of  attack),  and  not  the  details  of  the  motion  history  up  to  that  point  in 
time.  Using  this  idea,  the  knowledge  of  the  variation  of  separation  angle 
versus  pitch  rate  O*  for  constant  pitch  rate  motions  is  sufficient  to  construct  a  model 
to  predict  the  separation  angle  for  the  three-segment  trajectory  experiments.  The 
model  performs  remarkably  well  for  high  pitch  rates;  the  agreement  with  the 
experiments  is  somewhat  less  satisfactory  at  lower  range  of  pitch  rates. 

■  Shaping  the  pitch  trajectory  can  be  used  to  delay  the  onset  of  leading  edge 
separation  beyond  that  achieved  with  constant  pitch  rate  motions. 
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3.  Conclusions 


The  basic  physical  understanding  of  the  mechanisms  for  2-D  incipient  leading-edge 
separation  is  still  not  complete.  This  is  true  even  for  relatively  low  Reynolds  numbers  of 
order  10^  where  computational  results  are  available.  Due  to  lack  of  appropriate 
experimental  results,  the  adequacy  of  these  computations  is  hard  to  assess.  Most  needed 
are  the  details  of  the  flow  behavior  (e.g.  the  velocity  field)  within  the  boundary  layer  near 
the  leading  edge  and  the  evolution  of  the  reverse-flow  regions.  Such  boundary-layer 
resolved  measurements  are  outside  the  reach  of  traditional  measurement  techniques.  Novel 
diagnostics  based  on  molecular  tagging  techniques  offer  new  possibilities  for  these  types 
of  measurements  covering  chord  Reynolds  numbers  from  10^  to  10^. 

Another  issue  that  requires  more  attention  is  the  flow  behavior  for  variable  pitch  rate 
trajectories.  The  majority  of  studies  to  date  have  concentrated  on  either  sinusoidal  or 
constant  pitch  rate  motions.  Our  preliminary  results  indicate  that  the  onset  of  leading  edge 
separation  can  be  delayed  further  when  variable  pitch  rate  motions  are  considered.  The 
deliberate  shaping  of  the  pitch  trajectory  for  the  purpose  of  optimization  of  separation  delay 
is  one  way  to  manage  the  flow  and  aerodynamic  behavior  of  an  airfoil. 
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Figure  1.  Tangential  velocity  flow  reversal  locations  just  above  the  airfoil  surface  and  the 
vorticity  magnitude  contours  at  selected  angles  of  attack  (f2*  =  0.1).  Curve  labeled  (1)  and 
the  circle  symbol  indicate  the  reverse  flow  front  originating  at  the  trailing  edge.  Curve 
labeled  (2)  and  the  square  symbol  indicate  a  new  region  of  forward  flow. 
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Figure  2.  Tangential  velocity  flow  reversal  locations  just  above  the  airfoil  surface  and  the 
vorticity  magnitude  contours  at  selected  angles  of  attack  (Q*  =  0.4).  Curve  labeled  (1)  and 
the  circle  symbol  indicate  the  reverse  flow  front  originating  at  the  trailing  edge.  Curve 
labeled  (3)  and  the  triangle  symbol  indicate  the  appearance  of  a  new  independent  reverse 
flow  region.  Curve  (2)  and  the  square  symbol  indicate  the  region  of  forward  flow. 
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Initial  Acceleration  Effects  on  Flow 
Evolution  Around  Airfoils  Pitching 
to  High  Angles  of  Attack 

Manoochehr  M.  Koochesfahani*  and 
Vanco  Smiljanovskit 
Michigan  State  University, 

East  Lansing,  Michigan  48824 

Introduction 

IRFOILS  pitching  rapidly  to  high  angles  of  attack,  and 
the  accompanying  phenomenon  of  dynamic  stall,  were 
first  investigated  for  the  case  of  sinusoidal  motions.  The  pro¬ 
gress  in  this  area  is  reviewed  by  McCroskey'  and  Carr.^  More 
recently,  constant  pitch  rate  (i.e.,  ramp)  motions  have  been 
receiving  a  great  deal  of  attention  due  to  their  applicability  to 
supennaneuverable  aircraft.  Some  of  the  earlier  measure¬ 
ments  of  integrated  load’*’  have  now  been  complemented  by 
computational  results*  ’  and  high-resolution  surface  pressure 
measurements';  whole-field  velocity  information  is  also  be¬ 
coming  available.’ 

Investigations  into  airfoils  pitching  at  constant  rate  have 
typically  considered  airfoils  pitching  from  a  zero  angle  of 
attack  to  some  large  angle  a„„  well  beyond  the  static  stall 
angle.  It  has  been  established  that  a  major  parameter  govern¬ 
ing  the  now  behavior  is  the  nondimensional  pitch  rate 
n*  =  dC/(2l/.),  where  dt  is  the  constant  pitch  rate,  C  the 
chord,  and  f/.  the  freestream  speed.  It  is  clear  that,  in  reality, 
the  actual  motion  of  the  airfoil  must  deviate  from  the  ideal 
ramp  due  to  the  Hnite  acceleration  and  deceleration  periods 
imposed  by  the  damping  of  the  drive  system  and  the  response 
charaaeristics  of  the  airfoil.  The  computations  of  "nomi¬ 
nally”  constant  pitch  rate  motions  also  include  a  brief  initial 
acceleration  period.  To  our  knowledge,  a  systematic  investiga¬ 
tion  of  the  effects  of  initial  acceleration  on  the  flow  charaaer¬ 
istics  of  an  airfoil  pitching  to  high  angles  of  attack  has  not 
been  undertaken.  We  note  that  studying  such  effeas  can  pro¬ 
vide  not  only  further  insight  into  the  processes  of  vorticity 
generation  and  accumulation  on  unsteady  surfaces,  but  also 
clues  as  to  how  these  processes  may  be  modified  or  controlled 
by  the  deliberate  shaping  of  the  pitch-motion  trajeaory. 

In  the  present  experiments,  flow  visualization  is  used  to 
monitor  the  ons«  of  leading-edge  separation  and  the  subse¬ 
quent  dynamic  stall  vortex  development  as  the  initial  accelera¬ 
tion  is  systematically  varied  in  magnitude  and  duration 
through  the  acceleration  phase  to  constant  pitch  rate.  The 
work  presented  here  considers  the  case  of  incompressible  flow 
at  low-chord  Reynolds  numbers  and  relatively  high  pitch  rate. 
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b)  Example  of  actual  trajectory  executed  by  the  airfoil;  (1*  =0.4 

Fig.  1  Constant  pitch  rale  motion  with  finite  accdcmtion  and  decel¬ 
eration. 


Experimental  Setup 

The  experiments  were  performed  in  a  water  channel  (Engi¬ 
neering  Laboratory  Design,  Inc.)  with  a  test  seaion  60  x  60 
cm  in  cross  section  and  240  cm  in  length.  The  airfoil  was  an 
NACA  0012  with  a  uniform  chord  of  C  =  8  cm.  a  span  of 
b  =  45  cm.  and  was  fitted  on  both  ends  with  false  walls 
parallel  to  the  water  channel  walls.  For  the  results  described 
here,  the  freestream  seed  was  set  to  (/»  =  10  cm/s  resulting  in 
a  chord  Reynolds  number  of  8  x  lO’.  A  dc  servo  motor  and  a 
digitial  servo  controller  (Galil,  DMC-610)  were  used  to  pitch 
the  airfoil  about  the  quarter  chord.  The  pitch-motion  trajec¬ 
tory  staned  at  zero  angle  of  attack,  reached  the  desired  con¬ 
stant  pitch  rate  of  d  after  a  period  T,  of  constant  acceleration, 
and  stopped  at  the  final  angle  of  attack  of  60  deg  after  a 
period  of  constant  deceleration  (see  Fig.  la).  The  accelera¬ 
tion  and  deceleration  periods  were  kept  equal  in  this  work.  We 
characterize  the  pitch  trajectory  by  the  nondimensional  pitch 
rate  0*  and  an  acceleration  parameter  e  =  Ta/Tc,  where  is 
the  "ideal”  constant  pitch  rate  time  scale  needed  for  the 
motion.  Note  that  the  parameter  e  gives  an  indication  of  the 
fraaion  of  the  motion  time  used  for  acceleration/decelera¬ 
tion.  The  panicular  case  i  considered  in  this  study  correspond 
to  (D*  =  0.4;  e  =  0.6,  0.15,  0.037)  and  (Q*  =  0.2;  c  =  0.15, 
0.037).  Typical  examples  of  the  actual  pitch  trajectory  exe¬ 
cuted  by  the  airfoil  are  shown  in  Fig.  lb. 

The  evolution  of  the  flow  was  monitored  using  the  hydro- 
gen-bubble  technique  and  laser  sheet  illumination  at  the  airfoil 
midspan  location.  The  hydrogen-bubble  wire  was  placed  ap¬ 
proximately  I  mm  upstream  of  the  airfoil  leading  edge  and 
was  pulsed  at  20  Hz.  Flow  images  were  sensed  by  a  charge 
coupled  device  camera  at  a  rate  of  60  fields/s  with  an  exposure 
lime  of  2  ms/field  and  acquired  by  a  digital  image  acquisition 
system  (Recognition  Concepts,  Inc.,  TRAPlX-5500)  onto  a 
hard  disk  in  real  time. 
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Results  and  Discussion 

The  evolution  of  the  flowfieid  for  the  case  (fl*  =  0.4; 
e  =  0  15)  is  illustrated  in  Fig.  2.  For  each  picture,  the  instanta¬ 
neous  angle  of  attack  and  the  elapsed  time  from  the  start  of 
the  motion  are  indicated.  The  flow  follows  the  contour  of  the 
airfoil,  and  in  that  sense  appears  fully  attached,  up  to  an  angle 
of  attack  a  =  27  deg.  By  the  time  the  airfoil  has  reached 
a  =  34  deg,  leading-edge  separation  and  vortex  formation 
have  already  begun  with  a  visual  signature  of  a  perturbation 
near  the  leading  edge.  The  first  visual  indication  of  onset  of 
leading-edge  separation  and  vortex  formation  was  found  to 
occur  at  (7,  =  0.64  s;  a,  =  32  deg)  after  close  inspection  of  the 
image  time  sequence;  T,  and  a,  refer  to  the  elapsed  time  and 
angle  of  attack  at  the  onset  of  leading-edge  separation,  it 
should  be  noted  that,  based  on  the  image  acquisition  rate,  we 
estimate  the  value  of  a,  to  d:  1  deg  accuracy.  The  sequence  of 
pictures  in  Fig.  2  also  shows  how  the  leading-edge  vortex 
grows  in  time  and  evolves  into  the  dynamic  stall  vortex.  The 
flow  evolution  shown  in  this  figure,  including  the  formation 
of  multiple  large-scale  vortices  above  the  suction  surface,  is 
similar  to  the  known  computational  and  experimental  results. 

The  effect  of  initial  acceleration  period  on  the  onset  of 
leading-edge  separation  for  the  cases  considered  here  were 
determined  from  data  similar  to  Fig.  2.  The  results  are  sum¬ 
marized  in  Table  1. 

The  results  show  that  the  onset  of  leading-edge  separation  is 
delayed  to  a  larger  angle  of  attack  as  the  pitch  rate  increases. 
This  is  consistent  with  the  well-established  dynamic  stall  delay 
for  increasing  pitch  rate.''’  The  most  important  observation, 
however,  is  that  while  the  elasped  time  for  leading-edge  sepa¬ 
ration  and  vortex  formation  is  affected  by  initial  acceleration, 
the  angle  of  attack  where  this  occurs  remains  virtually  un¬ 
changed.  In  fact,  many  of  the  details  of  subsequent  flow 
development  are  also  nearly  the  same.  Figure  3  shows  a  com¬ 
parison  at  selected  angles  of  attack  for  two  cases  with  fl*  =  0.4 
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Fig.  2  Evolution  of  the  flowfield  on  the  airfoil  suction  surface 
(n'»0.4.  e-0.15). 


Fig.  3  Comparison  of  flow  evolution  for  two  different  initial  accel¬ 
eration  periods  |0‘  =0.4). 


Table  1  Time  and  angle  of  attack  at 
the  onset  of  leading-edge  separation 


7,.s 

a,,  deg 

C*  =  0.4;  e  =  0.6 

0.80 

31 

«•  =  0.4;  e  =  0.15 

0.64 

32 

!)•  =  0.4;  e  =  0.037 

0.57 

31 

n*  =  0.2:  e  =  0.15 

0.97 

24 

£)•  =  0.2;  e  =  0.037 

0.90 

25 

and  initial  accelerations  which  are  different  by  a  factor  of  4. 
Note  the  similarity  of  the  various  events  at  similar  angles  of 
attack;  there  is,  of  course,  a  time  shift  between  the  occurrence 
of  similar  events  in  the  two  cases,  as  indicated  earlier. 

The  downstream  convection  of  the  dynamic  stall  vortex  was 
determined  for  all  of  the  cases  from  image  sequences  similar  to 
Fig.  2.  The  vonex  downstream  position  X  (along  freestream 
direction)  was  estimated  using  the  center  of  the  nearly  circular 
region  interpreted  to  be  the  signature  of  the  dynamic  stall 
vonex  (e.g.,  see  Fig.  2).  The  results  are  referenced  to  the  fixed 
pitch  axis  location  Xp  and  are  plotted  in  Fig.  4.  Note  that  the 
differences  in  the  times  for  the  onset  of  leading-edge  separa¬ 
tion  among  the  various  cases  have  been  taken  into  account  by 
using  the  relative  time  (7-7,).  The  main  result  from  Fig.  4  is 
that  the  initial  acceleration  period  also  has  little  influence  on 
the  downstream  convection  of  the  dynamic  stall  vonex  for  the 
cases  studied.  Note  that  the  vonex  convection  speed 
(determined  from  a  least-squares  straight  line  fit  to  data)  is 
0.44  and  0.40  for  B*  =  0.4  and  0.20,  respectively. 

To  summarize,  the  results  presented  indicate  that  the 
elapsed  time  for  the  occurrence  of  various  events  is  affected  by 
the  initial  acceleration;  the  angle  of  attack  where  these  evei  s 
occur  is  practically  unchanged,  however.  This  includes  the 
onset  of  leading-edge  separation  and  many  of  the  details  of  the 
dvTiamic  stall  vonex  formation,  downstream  convection,  and 
its  interactions.  These  results  imply  that  a  convenient  accelera¬ 
tion  profile  can  be  selected  for  experimental  and  computa¬ 
tional  studies  without  seriously  impacting  the  dynamics  of  the 
unsteady  stall  process. 

.A.  scaling  argument  based  on  unsteady  inviscid  flow  results 
had  earlier  suggested that  if  the  initial  acceleration  time  scale 
is  sufficiently  shon  compared  to  the  flow  time  scale,  the  onset 
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Fit.  4  Oownsirtaa  convection  of  the  dynimk  stall  vonex.  Solid 
lines  represent  the  Ieast>s4|uares  line  fit  to  data. 


of  leading-edge  separation  might  be  expected  to  be  insensitive 
to  the  motion  history,  as  is  resported  here.  This  is  consistent 
with  the  constant  pitch  rate,  inviscid  flow,  lift-curve  slope 
results  in  Ref.  4  which  indicate  that  at  high  enough  pitch  rates 
motion  history  effects  will  saturate.  We  should  mention  that 
in  all  of  I  .e  cases  studied  here,  the  acceleration  period  was 
shon  enough  that  it  ended  before  the  onset  of  leading-edge 
separation.  Reference  10  further  suggests  that  initial  accelera¬ 
tion  effects  may  become  important  for  a  combination  of  low 
pitch  rate  Q*  and  long  acceleration  period  e. 

We  note  that  the  results  presented  here  only  address  the 
timing  of  the  various  events  in  the  flowfield  development.  We 
do  not  know  at  this  time  how  other  important  quantities  such 
as  the  surface  pressure  gradient,  surface  vonicity  flux,  and  the 
integrated  load  on  the  airfoil  are  affeaed  as  we  change  the 
acceleration  period.  These  questions  are  currently  being  ad¬ 
dressed  by  computing  the  flows  discussed  here  using  a  two-di¬ 
mensional  Navier-Stokes  solver.  Preliminary  computational 
results"  corroborate  the  present  findings. 

Acknowledgments 

This  work  was  supponed  by  the  Air  Force  Office  of  Scien¬ 
tific  Research  Grants  AFOSR-89-0417  and  AFOSR-89-0130. 
Charles  Gendrich  is  gratefully  acknowledged  for  his  help  with 
the  preparation  of  the  manuscript. 

References 

■McCroskey.  W.  J..  "Unsteady  Airfoils.”  Annua/  Review  of  Fluid 
Mechania.  Vol.  14.  1982.  pp.  283-311. 

^Can.  L.  W..  “Progress  in  Analysis  and  Prediaion  of  Dynamic 
Stall."  Journal  of  Airerafi,  Vol.  23.  No.  1.  1988.  pp.  6-l7. 

^Francis.  M.  S..  and  Keesce.  J.  E..  "Airfoil  Dynamic  Stall  Pnfor- 
mance  with  Large-AmpL’tude  Motions."  AlAA  Journal,  Vol.  23.  No. 
II.  1983.  pp.  1633-1639. 

'Jumper.  E.  J..  Shrcck.  S.  J..  and  Dimmick.  R.  L..  "Lift-Curve 
Characteristics  for  an  Airfoil  Pitching  at  Constant  Rate."  Journal  of 
Aircraft.  Vol.  24.  No.  10.  1987.  pp.  680-687. 

’Lorber.  P.  F..  and  Cana.  F.  O..  "Airfoil  Dynamic  Stall  at  Con¬ 
stant  Pilch  Rate  and  High  Reynolds  Number."  Journal  of  Aircraft, 
Vol.  23.  No.  6.  1988.  pp.  348-556. 

‘Visbai.  M.  R..  "On  the  Formation  and  Control  of  the  Dynamic 
Stall  Vonex  on  a  Pitching  Airfoil."  AlAA  Paper  91-0006.  Jan.  1991. 

'Chia.  K.  N..  Yang.  J..  Osswald.  C.  A.,  and  Chia.  U..  "Study  of 
the  Dynamic  Stall  Mechanism  Using  Simulation  of  Two-Dimensional 
Navier-Stokes  Equations."  AlAA  Paper  91-0346.  Jan.  1991. 

'Acharya,  M..  and  Metwally.  M.  H..  "Unsteady  Preuure  Field  and 
Vonicity  Production  over  a  Pitching  Airfoil.”  AlAA  Journal.  Vol. 
30.  No.  2.  1992.  pp.  403-411. 

'Shih.  C..  Lourenco.  L..  Van  Dommeien.  L..  and  Krothapalli.  A.. 
"Unsteady  Row  Past  an  Airfoil  Pitching  at  Constant  Rate."  AlAA 
Journal.  Vol.  30.  No.  3.  1992.  pp.  1 133-1 161. 

"’Koochesfahani.  M.  M..  Smiljanovski.  V.,  and  Brown.  T.  A., 


"Effect  of  Initial  Acceleration  on  the  Flow  Development  of  the  Row 
Field  of  an  Ah  foil  Pitching  at  Constant  Rale."  Proceedings  of 
S.ASA/AFOSR/ARO  Workshop  on  Physics  of  Forced  Unsteady  Sep¬ 
aration.  NASA  Ames  Research  Center.  April  17-19.  1990.  NASA  CP 
3144.  1992.  pp.  317-332. 

"Gendrich.  C.  P..  Koochesfahani.  M.  M..  and  Visbal.  M.  R.. 
"Initial  Acceleration  Effects  on  the  Flow  Field  Development  around 
Rapidly  Pitching  Airfoils."  AlAA  Paper  93-0438.  Jan.  1993. 


Decav  of  Aircraft  Vortices  near 

•> 

the  Ground 

Milton  E.  Teske*  and  Alan  J.  Bilanint 
Continuum  Dynamics.  Inc..  Princeton.  New  Jersey  08543 

and 

John  W.  Barryi 

U.S.  Department  of  Agriculture  Forest  Senice. 

Davis.  California  95616 

Introduction 

N  a  review  of  the  state  of  knowledge  of  aircraft  vonices. 
Donaldson  and  Bilanin'  include  a  discussion  on  the  effects  of 
atmospheric  turbulence  on  the  aging  of  vortex  pairs.  One  of  the 
models  suggested  may  be  w  ritten 


where  F  is  the  vonex  circulation  strength  as  a  fuiKtion  of  time  r.  To 
the  initial  vortex  circulation  strength  at  r  =  0.  b  the  decay  coeffi¬ 
cient.  g  the  ambient  turbulence  level,  and  s  the  aircraft  semispan. 

Donaldson  and  Bilanin’  invoke  simple  arguments  to  suggest 
that  the  decay  coefficient— out  of  ground  effea — may  take  a  value 
of  0.41.  When  deuils  about  the  wakes  of  large  airc^  were  first 
investigated  (in  the  early  1970sl.  more  sophisticated  models,  using 
various  techniques  such  as  second-order  closure  of  the  Reynolds 
stress  equations.*  were  developed  with  an  eye  toward  examining 
this  decay  behavior.  Recently,  for  a  completely  different  reason, 
the  problem  of  quantifying  the  decay  coefficient  has  produced  a 
sizable  field  study  and  data  examination. 

Discussion 

Between  1985  and  1991.  a  senes  of  aircraft  flybys  above  tower 
grids  instrumented  with  propeller  anemometers  pi^uced  a  large 
data  base  on  vortex  motion  near  the  ground.^  These  tests  were  per¬ 
formed  jointly  by  the  United  Sutes  Department  of  Agriculture 
Forest  Service  and  the  United  States  Army  at  several  sites  in  north¬ 
ern  California  (through  Project  WIND),  and  at  Dugway  Proving 
Ground.  Utah.  Anemometer  tower  grids  recorded  the  ambient  ver¬ 
tical  velocity  time  histories  as  various  aircraft  repeatedly  traversed 
normal  to  the  grid.  These  digitized  velocity  traces  produced  an  air¬ 
craft  wake  signature  that  could  be  used  to  infer  the  strength  and 
lateral  and  vertical  motion  of  the  aircraft  vortex  pairs  generating 
the  traces.  Table  I  summarizes  the  complete  data  set. 

The  location  and  strength  of  the  vortices  traversing  the  ane¬ 
mometer  tower  grid  may  be  inferred  by  a  least-squares  analysis 
defining  an  error  £  as 
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Computational  results  are  presented  to  show  how  the  acceleration  period  at  the  start  of  constant  pitch  rate  trajec¬ 
tories  alTects  the  dynamic  stall  process.  Large  amplitude  motions  were  studied  on  the  basis  of  an  NACA  0012  airfoil  pitch¬ 
ing  about  the  quarter-chord  axis.  Constant  (ON/OPF)  acceleration  profiles,  with  non-dimensional  acceleration  periods 
ranging  between  0.039  to  0.6,  and  large  pitch  rates  (0.1  ^  O*  a  dc/2U_  £  0.4)  were  used.  The  initial  acceleration  is 
observed  to  affect  the  integrated  loads,  surface  pressure  distribution,  and  tbe  evolution  of  reversed  flow  regions  Just  above 
the  surface  of  the  airfoil  only  during  the  acceleration  period  and  for  a  relatively  short  time  fx  s  tU_/c  ^  0.25)  afterwards. 
.\fter  that  time,  all  of  these  quantities  only  depend  on  the  instanuneous  angle  of  attack  for  a  given  pilch  rate.  These  con¬ 
clusions  are  consistent  with  previous  experimental  flow  visualization  observations.  The  onset  of  leading  edge  separation  at 
high  and  low  pitch  rates  is  shown  to  be  characterized  by  different  processes.  At  low  pitch  rates  leading  edge  separation 
occurs  after  (he  reversed  flow  front  originating  at  the  trailing  edge  has  reached  the  leading  edge.  At  higher  pitch  rates 
leading  edge  separation  and  the  upstream  progression  of  the  (trailing  edge)  reversed  flow  front  develop  independently. 


Nomenclature 

c  airfoil  chtx^d 

DSV  dynamic  stall  vonex 

DSVP  dynamic  stall  vonex  suction  peak  at  the  surface 

e  dimensionless  acceleration  rate.  T/T, 

LESP  leading  edge  sueden  peak 

M.  free  stream  Mach  number 

Re^  Reynolds  number  based  on  chord  length 

S  distance  along  the  airfoil  surface 

T,  acceleration  period 

Te  time  required  to  pitch  to  0®-60*  when  T,-»0 
U-  free  stream  velocity 

a  angle  of  attack 

o,^  a  at  the  onset  of  leading  edge  separation 
T  dimensionless  time.  tUjc 

Q*  dimensionless  pitch  rate.  (zc/2U_ 
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I.  Introduction 

The  dynamic  stall  process  has  been  under  inves- 
tigadon  for  over  a  decade  now,  and  significant  pro¬ 
gress  has  been  made  toward  understanding  the  physi¬ 
cal  processes  associated  with  rapidly  pitching  an  air¬ 
foil  beyond  its  stadc  stall  angle  of  attack.  The 
reviews  of  McCroskey  (1982),  Carr  (1988),  and  Visbal 
(1990)  provide  good  descripdons  of  the  dynamic  stall 
process  as  it  is  currently  understood.  As  the  dynamic 
stall  problem  was  first  recognized  for  its  relevance  to 
helicopter  performance,  the  first  modon  trajectories 
studied  were  sinusoidal  in  nature,  and  a  reduced  fre¬ 
quency  k  w  (ix:/2U.  was  defined  as  one  of  the  control¬ 
ling  parameters. 

More  recently,  constant  pitch  rate  modons  have 
received  a  great  deal  of  attendon  for  their  application 
to  supermaneuverable  aircraft.  These  trajectories  have 
an  advantage  over  sinusoidal  modons  for  those  study¬ 
ing  the  physics  of  the  problein,  since  they  offer  less 
complicated  experiments  with  which  to  study  the 
phenomena  associated  with  dynamic  stall.  For  con¬ 
stant  pitch  rate  motions  a  non-dimensional  pitch  rate. 
£2*  ■  dc/2U..  has  been  introduced  as  one  of  the 
governing  parameters.  It  is  clear  that  the  initial  part 
of  the  modon  trajectory  is  disdnguished  by  a  non¬ 
constant  pitch  rate  modon  due  to  the  finite  accelera- 
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lion  penod  which  is  characterisiic  ot  a  real  device. 
Very  little  is  known  about  the  influence  of  this  initial 
acceleration  period  on  the  subsequent  flow  develop¬ 
ment  around  airfoils  pitching  to  high  angles  of  attack. 

In  an  experimental  flow  visualization  study  by 
Koochesfahani  and  Smiljanovski  (1993),  the  initial 
acceleration  was  systematically  varied  in  order  to 
document  its  effects  on  flow  development  for  constant, 
relatively  high  pitch  rate  motions.  The  motion  trajec¬ 
tory  consisted  of  a  constant  acceleration  period  fol¬ 
lowed  by  a  constant  pitch  rate  segment  and  finally  a 
constant  deceleration  portion.  The  airfoil  pitched 
from  U’  to  60°;  see  Figure  1.  The  acceleration  effects 
were  described  in  tenns  of  a  non-dimensional 
acceleration  period,  e  a  T,/Tc,  where  T,  is  the 
acceleration  period  and  Tc  refers  to  the  time  to  pitch 
from  u°  to  60°  for  the  ideal  (i.e.,  T,  0)  constant 
pitch  rate  motion.  It  was  found  that  the  angle  of 
attack  at  the  onset  of  leading  edge  separation  was 
practically  independent  of  the  acceleration  period. 
Funheimore.  the  subsequent  flow  developments  after 
leading  edge  separation,  such  as  the  evolution  of  the 
dynamic  stall  vortex  and  its  convection  over  the  air¬ 
foil.  were  similarly  insensitive  to  the  details  of  the 
Stan  of  the  motion.  In  all  the  cases  studied.  leading 
edge  separation  occurred  after  the  initial  acceleration 
period  had  ended.  That  investigation,  which  was 
qualitative  and  entirely  based  on  flow  visualization, 
only  addressed  the  timing  of  the  various  events  in^the 
flow  field  development.  The  effects  of  initial 
acceleration  on  other  important  features  of  the  flow, 
such  as  the  integrated  loads  on  the  airfoil  and  the  sur¬ 
face  pressure  gradient,  have  not  yet  been  quantified. 

.Ml  computational  studies  of  nominally  constant 
pitch  rate  motions  also  contain  an  initial  acceleration 
phase.  However,  except  for  one  preliminary  investi¬ 
gation  (Visbal,  1986),  none  of  the  previous  studies  has 
specifically  addressed  the  effects  of  initial  acceleration 
on  the  dynamic  stall  process.  The  current  work, 
which  is  the  computational  counteipan  of  the  flow 
visualization  study  of  Koochesfahani  and  Smiljanovski 
(1993).  aims  to  quantify  the  initial  acceleration  effects 
on  the  pitching  airfoil  flow  field  using  two  dimen¬ 
sional  full  Navier-Stdces  computations.  These  effects 
are  characterized  here  in  terms  of  the  computed  flow 
visualizations,  the  integrated  load  on  the  airfoil,  the 
surface  pressure  distributions,  and  the  evolution  of 
reversed  flow  fronts  just  above  the  airfoil  surface. 

II.  Computational  Procedure 

The  current  computations  are  based  on  the 
Navier-Stokes  solver  of  Visbal  and  Shang  (1989). 


This  code  uses  the  implicit  approximately-factored 
algorithm  of  Beam  and  Warming  (1978)  to  solve  the 
unsteady  compressible  form  of  the  Navier-Stokes 
equations  written  in  the  strong  conservation  form. 
The  usual  forms  of  implicit  and  explicit  nonlinear  dis¬ 
sipation  as  described  by  Jameson.  Schmidt  and  Turkel 
(1981)  are  employed  in  order  to  ensure  numerical  sta¬ 
bility.  Further  details  can  be  found  in  Visbal  and 
Shang  (1989). 

The  computational  parameters  were  selected  to 
allow  comparison  with  the  range  of  experimental  con¬ 
ditions  in  the  water  tunnel  studies  of  Koochesfahani 
and  Smiljanovski  (1993)  and  Brown  (1992).  The 
computed  flow  field  was  that  of  an  NACA  0012  airfoil 
pitching  about  the  quarter-chord  axis  from  zero  to  60° 
angle  of  attack,  cl  The  freestream  Mach  number  was 
0.1,  and  the  chord  Reynolds  number  was  12.000.  The 
maximum  density  variation  was  less  than  3%  for  this 
Mach  number,  which  should  enable  acceptable  com¬ 
parisons  to  be  made  with  previous  water  tunnel  stu¬ 
dies.  Pitch  rates  in  the  range  0.1  <  £2*  <  0.4  were 
computed  with  acceleration  values  in  the  range 
0.001  <  e  <  0.6.  In  paper  results  will  be 
presented  only  for  £2*  =  '■  ..  0.2,  and  0.4  with  e  =  0.6, 
0.15,  and  0.039.  The  entire  computational  flow  field 
was  stored  every  Vi°  from  0*  through  50°.  The  flow 
was  expected  to  be  laminar,  at  least  during  the  initial 
phase  of  the  separation  process,  so  a  turbulence  model 
was  not  required. 

Particle  tracking  was  added  to  the  code  for  this 
study,  so  that  the  computational  results  could  be  com¬ 
pared  directly  with  experimental  hydrogen  bubble  flow 
visualizations.  A  fixed  number  of  "particles"  was 
introduced  into  the  computational  domain,  and  their 
positions  were  advanced  using  a  first-order  marching 
scheme.  The  particles  were  given  initial  locations  in 
space  such  that  they  simulated  lines  of  hydrogen  bub¬ 
bles  from  a  pulsed  wire.  Gendrich,  Koochesfahani. 
and  Visbal  (1992)  have  shown  a  very  good  correspon¬ 
dence  between  these  computed  "flow  visualizations" 
and  those  obtained  experimentally.  In  addition,  the 
visual  appearance  of  a  "bulge"  near  the  nose  of  the 
airfoil  at  the  onset  of  leading  edge  separation  corre¬ 
lates  strongly  with  changes  in  the  character  of  the  vor- 
ticity  field  and  the  tangential  velocity  just  above  the 
airfoil’s  surfac 


Results  and  Discussion 

Koochesfahani  and  Smiljanovski  (1993) 
observed  that  varying  the  initial  acceleration  period 
did  not  seem  to  affect  the  visual  appearance  of  the 
dynamic  stall  process.  In  Figure  2  we  present  a  com- 


pahson  between  their  experimental  and  our  computed 
tlow  visualizations  for  Q*  =  0.4  at  various  angles  of 
attack.  lUusuating  the  onset  of  leading  edge  separation 
*  for  two  different  accelerations.  There  are  many  more 
hydrogen  bubbles  in  the  experimental  flow  field  than 
particles  in  our  computed  domain,  which  accounts  for 
the  more  continuous  appearance  of  the  experimental 
visualizations.  Nevertheless,  the  computations  reveal 
the  presence  of  essentially  the  same  features  develop¬ 
ing  above  the  surface  of  the  airfoil.  At  IV  the  flow 
near  the  leading  edge  is  seen  to  follow  the  contour  of 
the  airfoil,  and  to  that  extent  exhibits  the  behavior  of 
an  attached  flow.  At  34®  a  "bulge"  has  formed  near 
the  leading  edge,  and  visually  we  could  identify  the 
onset  of  leading  edge  separation  to  have  occurred  at 
an  angle.  a,ep,  somewhere  between  27®  and  34®. 
Based  on  the  computed  flow  visualization  time  series, 
otjep  was  determined  to  be  approximately  33Vi®.  As 
the  airfoil  pitches  to  higher  angles  of  attack,  the  bulge 
develops  into  the  dynamic  stall  vonex  (DSV),  while 
secondary  vonex  rollup  and  pairing  occurs  over  the 
aft  ponion  of  the  airfoil.  The  computational  results 
shown  in  Figure  2  indicate  that  the  onset  of  leading 
edge  separation  is  not  influenced  by  the  initial 
acceleration,  in  agreement  with  the  experimental 
obser%'ations. 

To  quantify  the  initial  acceleration  effects,  we 
first  discuss  the  Iwhavior  of  the  integrated  loads  on  the 
airfoil.  The  computed  lift,  drag,  and  moment  curves 
for  fl*  =  0.4  are  shown  in  Figure  3.  The  spiky 
overshoots  near  the  beginning  of  each  curve  are 
believed  to  be  due  to  apparent  mass  effects.  The 
important  result  is  that  each  force  coefficient 
approaches  a  curve  which  depends  only  on  the  instan¬ 
taneous  angle  of  attack  (i.e..  is  independent  of  e) 
within  a  time  scale  of  t  =  0.26  after  ^e  end  of  the 
:icceleration  period  (see  Table  1);  x  is  the  non- 
dimensional  time  defined  by  x  s  tUJc.  For  the  par¬ 
ticular  pitch  rate  considered  in  Figure  3.  this  time 
scale  corresponds  to  a  maximum  change  in  the  angle 
of  attack  of  10®  beyond  the  end  of  acceleration.  The 
sudden  change  in  the  behavior  of  the  force  coefficients 
at  42®  for  the  case  of  e  =  0.6  is  due  to  the  beginning 
of  the  deceleration  period  for  that  case. 

Figures  4  and  S  show  the  lift,  drag,  and  moment 
coefficients  versus  angle  of  attack  for  £2*  =  0.2  and 
0.1.  respectively.  The  same  acceleration  effects  are 
observed  as  for  the  higher  pitch  rate,  namely,  witltin 
X  =  0.2S  after  the  end  of  the  acceleration  period,  the 
force  coefficients  have  become  independent  of 
xceleration  history.  Comparisons  with  the  experi¬ 
ments  of  Graham  and  Strickland  (1986).  Francis  and 
Keesee  (1985),  and  the  computations  of  Ghia.  Yang. 
Osswald,  and  Ghia  (1991)  are  included  in  Figures  4 


and  S.  Overall  agreement  with  the  experiments  is 
quite  good:  variations  are  probably  due  to  differences 
in  airfoil  type,  Reynolds  number,  and  iireestream  Mach 
number. 


e 

Acceleration 
ends  at: 

0.039 

0.15 

0.60 

a=  1.2® 
a  =  4.5® 
a  =  18.0® 

Table  1.  Angle  of  attack  at  which  acceleration  ends. 
Note  that  these  angles  are  independent  of  pi'  '^  rate 
for  constant  acceleration  motion  profiles.  De  on 
begins  at  42®  for  e  =  0.60. 

Since  the  integrated  loads  are  primarily  the 
result  of  the  surface  pressure  disuibution.  and  leading 
edge  separation  is  connected  to  the  surface  pressure 
gradient,  the  effects  of  xceleration  history  on  the 
development  of  the  surfxe  pressure  and  pressure  gra¬ 
dient  profiles  are  presented  next  The  typical  evolu¬ 
tion  of  these  profiles  at  a  fixed  value  of  initial 
acceleration,  e,  is  first  illustrated  in  Figures  6  and  7. 
Note  the  formation  of  the  leading  edge  suction  peak 
(LESP)  and  its  associated  adverse  pressure  gradient, 
followed  by  the  formation  of  the  dynamic  stall  vortex 
peak  (DSV?).^  The  first  indication  of  DSVP  formation 
occurs  between  16®  and  18.5*.  The  DSVP  magnitude 
is  initially  less  than  that  of  the  LESP,  but  it  grows  in 
strength  and  fitutlly  exceeds  the  leading  edge  suction 
peak.  While  the  LESP  remains  prxtically  fixed  at  a 
location  near  the  leading  edge,  the  DSVP  moves 
toward  the  trailing  edge  at  a  fairly  constant  rate.  The 
development  of  the  pressure  field  displayed  in  Figures 
6  and  7  agrees  very  well  with  that  presented  by 
Acharya  and  Metwally  (1992),  Visbal  (1991),  and  oth¬ 
ers. 

The  effects  of  varying  the  initial  xceleration  on 
the  development  of  the  surfxe  pressure  distribution 
are  illustrated  in  Figures  8-12.  Results  are  shown  for 
two  pitch  rates  at  selected  angles  before  and  after  the 
onset  of  leading  edge  separation.  The  evolution  for 
the  lower  pitch  rate  (£2*  =  0.1)  at  a  =  3®.  16®  (before 

and  20®  (after  o,,^  is  presented  in  Figures  8-10. 
At  3*  the  acceleration  has  already  ended  fenr 
e  =  0.039,  while  it  continues  for  the  other  two  e 
values.  The  pressure  over  prxtically  the  entire  surface 
is  different  for  exh  value  of  e.  with  the  exception  of 


'  The  DSVP  refers  to  the  pressuie  fninimum  at  the  airfoil 
surface,  and  not  the  pressure  mmimiim  in  the  core  of  the 
dynamic  stall  vonex.  This  dehnihoo  of  the  DSVP  is  the  same 
as  that  used  by  Adiaiya  and  Metwally  (1992). 
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the  stagnation  region.  Note,  however,  that  the  pres¬ 
sure  gradient  shows  much  less  sensitivity  to  variations 
in  e.  Similar  behavior  of  the  surface  pressure  gradient 
is  exhibited  at  a  =  16®.  just  before  the  formation  of 
the  DSVP:  note  that  by  this  angle  of  attack  the  surface 
pressure  profile  has  relaxed  to  a  single  curve  nearly 
independent  of  the  variations  in  e.  This  lack  of  sensi¬ 
tivity  by  the  surface  pressure  gradient  to  the  initial 
acceleration  history  is  responsible  for  the  correspond¬ 
ing  insensitivity  of  the  leading  edge  separation  angle 
observed  experimentally  and  computationally. 

After  acceleration  has  ended  for  all  three  cases 
and  leading  edge  separation  has  begun  (Figure  10. 
a  =  20®).  even  the  surface  pressure  field  exhibits  an 
insensitivity  to  e,  similar  to  that  shown  by  the  surface 
pressure  gradient.  This  similarity  among  the  different 
acceleration  cases  for  both  the  surface  pressure  and 
gradient  is  consistent  with  the  previous  experimental 
observations  indicating  that  the  subsequent  develop¬ 
ment  of  the  dynamic  stall  flow  field  is  not  influenced 
by  the  initial  acceleration  history  for  the  range  of 
parameters  investigated.  Figures  11  and  12  showing 
the  pressure  field  before  and  after  for  the  higher 
pitch  rate  of  Q*  =  0.4  suppon  all  the  assenions  made 
earlier.  For  this  pitch  rate  the  pressure  field  exhibits 
even  less  dependence  on  initial  motion  history  as 
compared  to  the  £2*  =  0.1  case.  It  should  be  noted 
that  the  discussion  presented  here  refers  to  the 
behavior  of  the  surface  pressure  field  over  the  first 
30%  of  the  chord  (see  Figs.  8-12).  There  are 
differences  in  the  pressure  field  in  the  aft  region  of  the 
airfoil  which  are  responsible  for  the  variation  in  the 
loads  among  the  different  cases;  see  Figs.  4  and  5. 
These  differences  do  not  have  a  significant  influence 
on  leading  edge  separation. 

The  discussion  above  addresses  the  flow 
behavior  for  only  selected  angles  of  attack.  The  evo¬ 
lution  of  the  surface  pressure  field  for  the  entire 
motion  has  been  analyzed  in  Gendiich,  Kooches- 
fahani.  and  Visbal  (1993)  in  terms  of  the  magnitude 
and  location  of  the  various  surface  pressure  minima 
and  pressure  gradient  maxima  corresponding  to  the 
LESP  and  DSVP  evolution.  The  main  conclusion 
from  these  results  (not  shown  here)  is  that  the  effects 
of  initial  acceleration  are  limited  to  a  short  period  (no 
more  than  t  =  0.2S)  after  the  end  of  acceleration.  In 
general,  the  development  of  the  pressure  gradient 
maxima  is  affected  less  than  the  evolution  of  the  pres¬ 
sure  minima,  and  the  minimum  pressure  magnitude  is 
less  sensitive  to  acceleration  history  than  its  location. 
After  a  certain  angle  of  attack  (a  =  6®  for  £2*  =  0.1. 
a  =  8®  for  £2*  =  0.2.  and  a  =  10®  for  £2*  =  0.4)  the 
magnitudes  and  locations  for  both  the  pressure  miruma 
and  pressure  gradient  maxima  become  independent  of 


initial  motion  history  and  depend,  for  each  pitch  rate, 
only  on  the  instantaneous  angle  of  attack. 

One  of  the  prominent  features  of  the  unsteady 
separation  process  on  a  pitching  airfoil  is  the  existence 
of  a  thin  region  of  reversed  flow  over  the  suction  sur¬ 
face  of  the  airfoil.  The  evolution  of  this  reversed  flow 
is  believed  to  be  intimately  connected  with  the  leading 
edge  separation  process.  The  effea  of  the  initial 
acceleration  on  the  development  of  various  regions  of 
reversed  tangential  flow  in  the  rotating  frame  of  the 
airfoil  is  characterized  in  Figure  13  for  three  different 
pitch  rates.  For  each  case,  the  progression  of  the 
reversed  flow  originating  at  the  trailing  edge  has  been 
labeled  (1).  and  that  which  results  from  the  voracity 
concentration  corresponding  to  the  dynamic  stall  vor¬ 
tex  has  been  labeled  (2).  Other  regions  of  flow  rever¬ 
sal  are  also  shown  in  this  figure;  they  correspond  to 
other  concentrations  of  vonicity  (vonex  roll-ups) 
existing  over  the  upper  surface. 

Figure  13  (a).  £2*  =  0.1.  illustrates  that  the  front 
of  the  reversed  flow  region  labeled  (1)  begins  near  the 
trailing  edge  and  moves  upstream.  When  this  front 
nears  the  leading  edge  (a  =  20®).  a  new  region  of  for¬ 
ward  flow  labeled  (2)  forms,  marking  the  occurrence 
of  leading  edge  separation.  It  is  very  important  to 
note  that  this  behavior  is  quite  different  at  higher 
pitch  rates;  sec  for  example  Figure  13  (c)  for 
£2*  =  0.4.  In  this  case  the  progression  of  the  reversed 
flow  front  originating  at  the  nailing  edge  never  nears 
the  leading  edge  (it  only  reaches  S/c  =  0.3)  before 
leading  edge  separation  is  observed  (as  indicated  by 
contour  (2)).  For  this  higher  pitch  rate  a  separate 
region  of  reversed  flow  labeled  3  develops  near  the 
leading  edge  at  approximately  20®.  Experimental  evi¬ 
dence  for  this  occurrence  of  leading-edge  reverse  flow 
independent  of  the  trailing-edge  reverse  flow  has 
recently  been  provided  by  Acharya  (private  communi¬ 
cation). 

As  far  as  the  initial  acceleration  effects  are  con¬ 
cerned,  a  longer  acceleration  period  (e.g.,  e  =  0.6) 
appears  to  speed  the  growth  of  the  reversed-flow 
region  from  the  trailing  edge  toward  the  leading  edge. 
However,  the  important  result  from  figure  13  is  that 
regardless  of  which  of  the  two  processes  described 
above  is  characteristic  of  the  flow  (i.e..  leading  edge 
stall  versus  trailing  edge  stall),  the  initial  acceleration 
does  not  significantly  influence  that  process. 


IV.  Conclusions 

The  unsteady  flow  field  around  an  NACA  0012 
airfoil  pitching  to  large  angles  of  attack  was  presented 
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for  three  different  pitch  rates,  fl*  =  0.1.  0.2.  and  0.4. 
;tnd  three  different  initial  constant  acceleration 
periods,  e  =  0.039.  0.15.  and  0.6.  These  are 
moderately  large  pitch  rates  for  which  leading  edge 
separation  begins  after  the  acceleration  has  ended:  no 
cases  have  been  considered  for  which  the  airfoil  was 
still  accelerating  at  the  onset  of  leading  edge  separa¬ 
tion.  The  contputational  parameters  were  selected  to 
allow  comparison  with  earlier  experimental  studies. 

The  computed  flow  visualizations  support  the 
conclusions  reached  by  the  experimental  obsmations 
of  Koochesfahani  and  Smiljanovski  (1993)  that  the 
onset  of  leading  edge  separation  and  the  subsequent 
development  of  the  dynamic  stall  vortex  is  not 
influenced  by  the  initial  motion  history  for  the  range 
of  parameters  studied.  The  present  computations  pro¬ 
vide  detailed  quantitative  information,  not  available 
from  those  experiments,  about  the  integrated  loads, 
surface  pressure  distributions,  and  the  evolution  of 
reversed  flow  regions  just  above  the  surface  of  the  air¬ 
foil.  The  results  show  that  the  initial  acceleration 
impacts  these  quantities  only  during  the  acceleration 
period  and  for  a  relatively  shon  time  (x  5  0.25)  after¬ 
wards.  After  that  time,  all  of  these  quantities  only 
depend  on  the  instantaneous  angle  of  attack  for  a 
given  pitch  rate.  These  findings  are  in  qualitative 
agreement  with  the  inviscid  flow  lift-curve  slope 
results  of  Jumper,  Shreck.  and  Dimmick  (1987)  which 
indicate  that  at  high  enough  pitch  rates  motion  history 
effects  will  saturate.  We  draw  attention  to  the  impor¬ 
tant  implication  of  our  results  that  a  convenient 
acceleration  profile  can  be  selected  for  experimental 
and  computational  studies  without  seriously  impacting 
the  dynamics  of  the  unsteady  stall  process. 

The  results  presented  also  demonstrate  that  two 
different  processes  characterize  the  onset  of  leading 
edge  separation  at  high  and  low  pitch  rates.  At  low 
pitch  rates  the  onset  of  leading  edge  separation  occurs 
after  the  reversed  flow  firont  originating  at  the  trailing 
edge  has  reached  the  leading  edge.  At  higher  pitch 
rates,  these  two  features  develop  independently;  lead¬ 
ing  edge  separation  occurs  without  the  reversed  flow 
front  from  the  trailing  edge  having  ever  reached  the 
leading  edge  region. 
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Figure  1.  Ideal  and  typical  motion  profiles.  esT^/Tf.  and  Q*  a  dc/2U, 


Figure  2.  Experimental  and  computational  flow  visualizations  for  0.4  and  two  different  acceler 
adons.  a)  experiments,  e=0. IS;  b)  computations,  e=0. IS;  c)  experiments,  e=0.60;  andd) 
computations,  e=0.60.  Experimental  flow  visualizations  are  from  Koochesfahani  and  Smil- 
janovski  (1993). 
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Figure  4.  a)  Lift,  b)  Drag,  and  c)  Moment  coefficients.  Q*=  0.2.  . e  =  0.039. 

0.15. - e  =  0.60.  and  ♦ - •  Graham  and  Strickland  (1986) 
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Figure  6.  Development  of  the  surface  pressure  profile  during  a  dynamic  stall  event,  Q*  =  0. 1,  e  = 
0.15;  a)  inception  and  b)  growth  of  the  dynamic  stall  vortex  peak  (DSVP) 
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Figure  7.  Surface  pressure  gradient  evolution  during  a  dynamic  stall  event.  Q*  =  0. 1,  e  =  0.15. 
Angles  are  the  same  as  in  Figure  6  a)  and  b). 
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Figures,  a)  Surface  pressure  and  b)  surface  pressure  gradient  profiles  at  3°.  Q*=0.1.  Acceler¬ 
ation  has  ended  for  only  the  e  =  0.039  trajectory.  . e  =  0.039.  _ e  =  0.15.  and 
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Figure  9.  a)  Surface  pressure  and  b)  surface  pressure  gradient  profiles  just  before  the  onset  of 
leading  edge  separation  (16°").  =0.1.  . e  =  0.039.  - e  =  0.15.  and - « 


Figure  10.  a)  Surface  pressure  and  b)  surface  pressure  gradient  profiles  just  after  the  onset  of 
leading  edge  separation  (20°),  Q*  =0.1.  . e  =  6.039. - e  =  0.1 5.  and - e  = 
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Figure  11.  a)  Surface  pressure  and  b)  surface  pressure  gradient  profiles  just  before  the  onset  of 
leading  edge  separation  (22'^).  Q*  =  U.4 . e  =  0.039,  - e  =  0. 15,  and - e 
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Figure  12.  a)  Surface  pressure  and  b)  surface  pressure  gradient  profiles  just  after  the  onset  of 

leading  edge  separation  ( 30*^),  Q*  =  0.4 . e  =  0.039,  - e  =  0. 1 5,  and - e  = 

0.60.  " 
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Figure  13.  Tangential  velocity  flow  reversal  locations  just  above  the  airfoil  surface;  a)  =  0. 1 
b)  n*  =  U.2;  andc)  Q*=U.4.  o  e  =  0.039,  oe  =  0.15,  and  ve=0.60.  See  text  for  a 
description  of  the  curves  marked  1, 2,  and  3. 
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